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Abstract We performed fully self-consistent calculations of p-nuclear bound 
states within the relativistic mean-field (RMF) model. The G-parity moti¬ 
vated p-meson coupling constants were adjusted to yield potentials consistent 
with p-atom data. We confirmed large polarization effects of the nuclear core 
caused by the presence of the antiproton. The p absorption in the nucleus 
was incorporated by means of the imaginary part of a phenomenological opti¬ 
cal potential. The phase space reduction for the p annihilation products was 
taken into account. The corresponding p width in the medium significantly 
decreases, however, it still remains considerable for the p potential consistent 
with experimental data. 
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1 Introduction 

The study of the antiproton-nucleus interaction is an interesting issue which 
has attracted renewed interest in recent years at the prospect of future exper¬ 
iments at the FAIR facility. The p-nuclear bound states and the possibility 
of their formation have been studied in refs. m- These considerations are 
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supported by a strongly attractive potential that the p feels in the nuclear 
medium. Within the RMF approach the real part of the p-nucleus poten¬ 
tial derived using the G-parity transformation is ReV^ ~ —650 MeV deep at 
normal nuclear density. However, the experiments with p atoms [6] and p scat¬ 
tering off nuclei at low energies FT] favor shallower real part of the p-nucleus 
potential in the range of —(100 — 300) MeV in the nuclear interior. An im¬ 
portant aspect of the p-nucleus interaction is p annihilation which appears 
to be the dominant part of the interaction. Nevertheless, the phase space for 
the annihilation products should be significantly suppressed for the antiproton 
bound deeply in the nuclear medium, which could lead to the relatively long 
living p inside the nucleus pH. 

In this contribution, we report on our recent fully self-consistent calcu¬ 
lations of p-nuclear bound states including p absorption in a nucleus. The 
calculations are performed within the RMF model [3]. Dynamical effects in 
the nuclear core caused by the antiproton and the phase space suppression for 
the p annihilation products are studied for various nuclei. 

In Section 2, we briefly introduce the underlying model. Few selected rep¬ 
resentative results of our calculations are discussed in Section 3. 


2 Model 


The p-nucleus interaction is described within the RMF approach. The inter¬ 
action among (anti)nucleons is mediated by the exchange of the scalar (cr) and 
vector (w M , p p ) meson fields, and the massless photon field A The standard 
Lagrangian density Cn for nucleonic sector is extended by the Lagrangian 
density C P describing the antiproton interaction with the nuclear medium (see 
ref. [J] for details). The variational principle yields the equations of motion for 
the hadron fields involved. The Dirac equations for nucleons and antiproton 
read: 

[-iaV + 0(m j + S j ) + V j ]i/>f = e^f, j = N,p , (1) 


where 


Sj — (Jrrjj ( 7 , 


Vj = g uj uj 0 + g P jPoT3 + ej 


1 + t 3 
2 


Aq 


( 2 ) 


are the scalar and vector potentials. Here, rrij stands for (anti)nucleon mass; 
g<rj , guij , g P j, and ej are the (anti)nucleon couplings to corresponding fields, 
and a denotes single particle states. The Klein-Gordon equations for the boson 
fields acquire additional source terms due to the presence of p: 


(-A + m 2 a + g 2 cr + g 3 cr 2 )cr = -g a NPSN - g^pPsp , 
(-A + m 2 u + dwo)w 0 = g^NpvN + g^pPvp , 
(—A + m 2 p )po = g pN P in + g P pPip , 


—AA 0 = ej sspqn + epPQp , 


( 3 ) 


where psj, Pvj, Pij and pQj are the scalar, vector, isovector, and charge den¬ 
sities, respectively, and rn a . rri ^, rn p are the masses of the considered mesons. 
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In this work, the nucleon-meson coupling constants and meson masses were 
adopted from the nonlinear RMF model TM1(2) [5] for heavy (light) nuclei. 
The system of the coupled Dirac CD and Klein -Gordon j3j equations is solved 
fully self-consistently by iterative procedure. 

In the RMF model, the nucleon in a nucleus moves in mean fields created 
by all nucleons, i. e., the nucleon feels repulsion as well as attraction also from 
itself. In ordinary nuclei this nucleon self-interaction has only a minor (1/A) 
effect. However, the potential acting on p in a nucleus is much deeper and the 
impact of the p self-interaction could become pronounced. In order to exclude 
this unphysical p self-interaction we omitted the antiproton source terms in 
the Klein-Gordon equations for the boson fields acting on the 

The p-nucleus interaction is constructed from the p-nucleus interaction 
with the help of the G-parity transformation: the vector potential generated 
by the w meson exchange thus changes its sign and becomes attractive. As a 
consequence, the total p potential will be strongly attractive. However, the G- 
parity transformation should be regarded as a mere starting point to determine 
the p-meson coupling constants. Various many-body effects, as well as the 
presence of strong annihilation channels could cause significant deviations from 
the G-parity values in the nuclear medium. Therefore, we introduce a scaling 
factor £ S (0,1) for the p-meson coupling constants [T|: 

9c rp = £ 9<?n j 9cjP = £ 9cjN j 9pp = £ 9pN ■ (4) 


The p annihilation in the nuclear medium is described by the imaginary 
part of the optical potential in a ‘tp’ form adopted from optical model phe¬ 
nomenology {BJ: 

2plmV opt (r) = -47t (l + J!— A I mb 0 p(r) , (5) 

v m N A J 


where p is the p-nucleus reduced mass. While the density p(r) was treated 
as a dynamical quantity evaluated within the RMF model, the parameter 
Im6o = 1.9 fm was determined by fitting the p atom data [6]. 

The energy available for the p annihilation in the nuclear medium is usually 
expressed as y/s = mp+mN — Bp — Bn, where Bp and Bn is the p and nucleon 
binding energy, respectively. The phase space available for the annihilation 
products is thus considerably suppressed for the deeply bound antiproton. 

The phase space suppression factor f s for two body decay is given by [8] 


M 1 2 j [s-(wi+w 2 ) 2 ][s-(mi-m 2 ) 2 ] r- 

s y [M 2 — (to i + 771-2)] [Af 2 — (toi — TO.2) 2 ] " 


mi - 777-2) , 


( 6 ) 


where toi, m 2 are the masses of the annihilation products and M = to-p + to-at. 
For channels containing more than 2 particles in the final state the f s was 
evaluated with the help of the Monte Carlo simulation tool PLUTO [5]. 


1 It is to be noted that the self-interaction is directly subtracted in the Hartree—Fock 

formalism. 
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Fig. 1 The p density distribution in 208 Pb, 
calculated for different £ in the TM1 model 
with (left) and without (right) the p self¬ 
interaction. 



Fig. 2 The nuclear core density distribu¬ 
tion in 40 Ca and 40 Ca+p for £ = 0.2, calcu¬ 
lated in the TM1 model. The p density dis¬ 
tribution in 40 Ca+p (dotted line) is shown 
for comparison. 


3 Results 

The formalism introduced above was employed in the self-consistent calcula¬ 
tions of the p bound states in selected nuclei. First, we did not consider the p 
annihilation in the nuclear medium and focused on the study of the dynam¬ 
ical effects caused by the presence of the antiproton in the nucleus. During 
our calculations we noticed a pronounced effect of the p self-interaction on 
the calculated observables. In Fig. |I] we present the p density distribution in 
208 Pb, calculated dynamically in the TM1 model for different values of the 
scaling parameter £. The central p density pp{ 0) calculated including the p 
self-interaction (left panel) reaches its maximum for £ ss 0.5 and then starts 
to decrease. It is due to the interplay between the negative value of Sp — Vp 
(absolute value of which increases with £), the p single particle energy, and 
the p rest mass, which affects the solution of the Dirac equation for the p wave 
function. On the other hand, when the p self-interaction is subtracted (right 
panel), the scalar and vector potentials are of comparable depth and the pp{ 0) 
increases gradually with £ and saturates at much higher values. It is to be 
stressed that the effect of the p self-interaction is negligible for £ < 0.2, which 
includes the p potentials consistent with p atom data. 

Our calculations revealed large polarization of the nuclear core caused by 
the p in the nuclear Is state. The nuclear core density in a p nucleus reaches 
2 — 3 times the nuclear matter density as illustrated for 40 Ca in Fig. [2] The p 
is localized in the center of the nucleus and the resultant nuclear core density 
is substantially enhanced over a small region, r < 1.5 fm. 

The nucleon single particle energies are affected by the presence of the p as 
well. Consequently the total binding energies of p nuclei increase considerably, 
B = -344.3 MeV for 40 Ca and B = -485.4 MeV for 40 Ca+p (£ = 0.2). 
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Fig. 3 The phase space suppression fac¬ 
tors f s as a function of the c.m. energy 

VS- 



Fig. 4 Upper components g (left) and the 
relation 0 for lower components / (right) of 
the lpi /2 anf l 1 Ps /2 P wave functions in 16 0, 
calculated dynamically in the TM2 model us¬ 
ing the real (no abs) and complex (+ abs) p 
potential. 


In order to account for p annihilation, we performed calculations using 
the complex potential presented in Section 2. We considered the suppression 
of the phase space for the annihilation products. In Fig. [3j the phase space 
suppression factors for the annihilation channels involved are presented as a 
function of yfs. As yfs decreases due to the p and N binding energies many 
channels become strongly suppressed or even closed. Moreover, the p-nucleon 
annihilation takes place in the nuclear medium. Therefore, the momentum 
dependent term in the Mandelstam variable s = (En + Ep) 2 — (pjy + Pp) 2 , 
where Ej = rrij — Bj , is non-negligible in contrast to the two body frame 
m- Our self-consistent evaluation of yfs including pp and pvr leads to an 
additional downward energy shift overlooked by many previous calculations. 

In Table [TJ we present the Is p single particle energies and widths in 16 0+p, 
calculated using the real and complex potentials consistent with p-atom data 
(£ = 0.2). To illustrate the role of the suppression factors f s we show the 
results of calculations without f s (‘Complex’), as well as including f s for yfs 
due to Bp and Bn (‘Complex+/ s ’) and for yfs with the additional downward 
energy shift due to the momenta of annihilating partners (‘+ys(p)’). The 
static calculations, which do not account for the core polarization effects, give 
approximately the same values of the p single particle energy for all cases. 
The single particle energies calculated dynamically are larger, which indicates 
that the polarization of the core nucleus is significant. When the phase space 
suppression is taken into account the p width is reduced by more than twice 
(compare ‘Complex’ and ‘Complex+/ s ’ in the last row of Table [l]). When 
treating yfs self-consistently including the p and N momenta, the p width is 
reduced by additional ~ 50 MeV, but still remains sizeable. The corresponding 
lifetime of the p in the nucleus is ~ 1 fm. 
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Table 1 The Is single particle energies ep and widths Tp (in MeV) in ie O+p, calculated 
dynamically (Dyn) and statically (Stat) within the TM2 model using the real and complex 
potentials consistent with p-atom data (see text for details). 



Real 

Complex 

Complex 

+ fs 

+V s (p) 


Dyn 

Stat 

Dyn 

Stat 

Dyn 

Stat 

Dyn 

Stat 

£p 

-193.7 

-137.1 

-175.6 

-134.6 

-190.2 

-136.1 

-191.6 

-136.3 

r P 

- 

- 

552.3 

293.3 

232.5 

165.0 

179.9 

144.7 


Finally, we discuss spin symmetry in p spectra. Static calculations with 
a real potential m revealed that antinucleon spectra in nuclei exhibit spin 
symmetry. We explored the p spectra considering the dynamical effects as 
well as the p absorption in the nucleus. In Fig. [4] the real parts of the upper 
components (left) and the differential relation jT2] 

+ = (J^~ fnr,e-i/ 2 (r) ( 7 ) 

for the real parts of the lower components (right) of the p wave function in 
lpi /2 and IP 3/2 states in 16 0 are plotted for £ = 0.2. We found that spin 
symmetry is well preserved in the p spectrum calculated fully self-consistently 
using the p potential consistent with the p atom data. The p annihilation 
causes only minor deviations. When different values of £ for the scalar and 
vector potentials are considered, spin symmetry holds only approximately. The 
deviations gradually increase with increasing difference between the scalar and 
vector potentials. 
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